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ABSTRACT 



A sample of large northern Spitzer Infrared Nearby Galaxies Survey (SINGS) galaxies was observed with the Westerbork Synthesis 
Radio Telescope (WSRT) at 1300 - 1760 MHz. In Paper II of this series, we described sensitive observations of the linearly polarized 
radio continuum emission in this WSRT-SINGS galaxy sample. Large-scale magnetic field structures of two basic types are found: 
(a) disk fields with a spiral topology in all detected targets; and (b) circumnuclear, bipolar outflow fields in a subset. Here we explore 
the systematic patterns of azimuthal modulation of both the Faraday depth and the polarized intensity and their variation with galaxy 
inclination. A self-consistent and fully general model for both the locations of net polarized emissivity at 1 - 2 GHz frequencies and 
the global magnetic field topology of nearby galaxies emerges. Net polarized emissivity is concentrated into two zones located above 
and below the galaxy mid-plane, with the back-side zone suffering substantial depolarization (by a factor of 4 - 5) relative to the 
front-side zone in its propagation through the turbulent mid-plane. The field topology which characterizes the thick-disk emission 
zone, is in all cases an axisymmetric spiral with a quadrupole dependance on height above the mid-plane. The front-side emission is 
affected by only mild dispersion (lO's of rad m"^) from the thermal plasma in the galaxy halo, while the back-side emission is affected 
by additional strong dispersion (lOO's of rad m"-) from an axi-symmetric spiral field in the galaxy mid-plane. The field topology in 
the upper halo of galaxies is a mixture of two distinct types: a simple extension of the axisymmetric spiral quadrupole field of the 
thick disk and a radially directed dipole field. The dipole component might be a manifestation of (1) a circumnuclear, bipolar outflow, 
(2) an in situ generated dipole field, or (3) evidence of a non-stationary global halo. 

Key words. ISM: magnetic fields - Galaxies: magnetic fields - Radio continuum: galaxies 



1. Introduction 

The magnetic fields in spiral galaxies are an important compo- 
nent, but their basic three dimensional topology remains largely 
unknown. Two of their main characteristics are however, known. 
First, the fields in relatively face-on spiral galaxies are seen to 
follow the spiral pattern traced in the optical morphology. In 
the handful of more edge-on galaxies that have been imaged 
to date, the field distributions are seen to extend into the halo 
regions, and have a characteristic X-shaped morphology (eg. 
iHeesen et al. 2009). Apart from these basic properties, the de- 
tails of the magnetic field topology are unknown. 

Observations of polarized flux, polarization vector orienta- 
tions, and Faraday rotation measures all provide information 
about the magnetic field associated with different electron pop- 
ulations and at different projections with respect to the line of 
sight. Synchrotron emission originates in ultrarelativistic elec- 
trons spiralling around magnetic field lines, is beamed in the di- 
rection of motion of the electron, and is polarized perpendicular 
to the orientation of the field line. Polarized synchrotron radi- 
ation and polarization vector orientation are thus direct tracers 
of the magnetic fields perpendicular to the line-of-sight (LOS), 
B±, within the region where both ordered magnetic fields and 
relativistic electrons are maximized. The Faraday rotation mea- 
sure (RM), or more generally the Faraday depth, O, that pertains 



to a given component of polarized emission, is sensitive to the 
integrated product of magnetic field component parallel to the 
LOS (B||) and the thermal electron density in the foreground of 
a polarized emission component: 
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The Faraday depth is defined to be positive when B points to- 
ward the observer, and negative when B points away. When 
assessing the magnetic field geometry traced by these observa- 
tional characteristics, it is essential to keep in mind that the ob- 
servable attributes may originate in distinct regions of space. The 
classical Faraday rotation measure (RM) is an observable quan- 
tity derived from the polarization angle difference(s) Ax between 
two (or more) frequency bands as RM - Ax/(A^ - A^). The em- 
pirically determined RM is only equivalent to the Faraday depth 
<S> for a simple background emitter plus foreground dispersive 
screen geometry. 

Polarized emission can become depolarized in a number of 
ways: beam depolarization can arise because the spatial resolu- 
tion element is large relative to the size of significant variations 
in the field orientation or the thermal electron content, while 
Faraday depolarization can arise because synchrotron emission 
and Faraday rotation take place in the same extended volume 
along the LOS. Polarized emission from different locations (ei- 
ther separated spatially or along the LOS) is affected by differ- 
ent amounts of Faraday rotation, such that at a given wavelength 
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there may be orthogonal polarization angles that cancel, yielding 
no net polarization at that wavelength. Beam depolarization can 
be circumvented in principle by using higher angular resolution, 
although the brightness sensitivity may then be insufficient to de- 
tect the extended emission at all. Faraday depolarization can be 
circumvented in principle by achieving a sufficiently complete 
sampling of the measurement space (relevant for measure- 
ments of RM and O), since cancellation effects are confined to 
discrete wavelengths or ranges of wavelength. 

All of these observables can be used to constrain the likely 
magne tic field topo logy in the galaxies observed. In a previous 
paper (' Heald et"al] (I2009D . hereafter Paper II), we presented our 
polarimetric results for a large sample of nearby galaxies ob- 
served to a comparable sensitivity limit of about 10 yidy beam"' 
RMS. In this paper, we begin by briefly summarizing the obser- 
vations and data reduction steps of Paper II in Sect. |2] Trends 
noted in the data are described in Sect. [3] In Sect. |4] we then 
explore how particular magnetic field geometries might relate to 
the observations. We conclude the paper in Sect.|5] 

2. Summary of observations and data reduction 

The observational parameters and data reduction techniques 
of the WSRT-SING S survey were presented in detail both by 
iBraun etaP (l2007l) . and specifically regarding the polarization 
data in Paper II. Here we recap the most important details. 
For more information, the reader is referred to Sect. 2 and 
Appendix A of Paper II. 

The data used in this analysis were obtained using the 
Westerbork Synthesis Radio Telescope (WSRT). Two observ- 
ing bands were used: of 1300 - 1432 MHz and 1631 - 1763 
MHz (centered on 22- and 18-cm, respectively), there being 
512 channels in each band and in all four polarization prod- 
ucts. Each galaxy in the WSRT-SINGS sample (refer to Paper 
II) was observed for 12 hr in the Maxi-short configuration of the 
WSRT. During each 12 hr synthesis, the observing frequency 
was switched between the two bands every 5 min. This provided 
an effective observation time of 6 hr per band, and good uv cov- 
erage in both bands. 

The data for both bands were analyzed using the 
Rotation Measure Synthe sis (RM-Synthesis) technique 
dBrentiens & de Bruvnll2005l see also Paper II). This provides 
the possible reconstruction of the intrinsic polarization vectors 
along each LOS, within the constraints set by the observing 
frequencies. The output of the RM-Synthesis procedure was 
deconvolved along the Faraday depth (O) axis, as described 
in Paper II. Polarized fluxes, polarization angles, and Faraday 
depths were extracted from these data and are discussed for 
each target galaxy in Paper II. In that paper, we also estimate the 
contribution to the RM from the Milky Way foreground using 
only background radio sources in the observed fields, rather 
than the target galaxies themselves. With this collection of data, 
we noted several patterns in the target galaxies, and that the 
basic patterns were common to the sample galaxies collectively. 
In this paper, we seek to explain these patterns using a common 
global magnetic field topology. 

3. Observational trends 

Several interesting patterns emerge from our study of a large 
sample of galaxy types. Polarized emission is found to origi- 
nate both in the disks of actively star-forming galaxies and in 
what appear to be either AGN or star-formation driven (circum- 
)nuclear or galactic wind outflows. "Disk" emission is relatively 



planar and detected out to large radii, whereas apparent "out- 
flow" components extend much further from the plane and are 
only detectable at small radii. The magnetic field orientations 
are in all cases simply related to the morphology of these com- 
ponents (as shown in Fig. 4 of Paper II). 

Disk fields have a spiral morphology that is strongly cor- 
related with the orientation and pitch angle of traditional trac- 
ers of spiral arms, such as massive stars and dust lanes. This 
is despite the polarized emission possibly being both coincident 
with, and independent of, other spiral arm tracers. Classic ex- 
amples of disk-dominated fields can be seen in galaxies such as 
NGC 628, 5194, and 6946, although they are present in almost 
all of our targets with varying detectability. On the other hand, 
outflow-related fields are typically oriented along the periphery 
of the bipolar lobes, and are often brightest close to the disk. The 
contribution of outflow-related components may be apparent in 
NGC 4569, 4631, and possibly 4736. 

Here, we illustrate and describe some of the trends present in 
the WSRT-SINGS dataset. We consider the polarized flux distri- 
bution (Sect. 13.1b . the Faraday depth distribution (Sect. 13.2b . the 
effects of depolarization in (Sect. [33] and Sect. 13.4b . In Sect.|4] 
we explore the predictions of a variety of global magnetic field 
topologies that can be used to model these trends. 

3.1. Polarized flux 

A remarkable pattern emerges for our sample relating to the ba- 
sic distribution of polarized intensity in galaxy disks at GHz fre- 
quencies. For small inclination angles, there is a general gradient 
in the average polarized intensity that is approximately aligned 
with the major axis of the target galaxy. This gradient from high 
to low polarized intensity has the same sign in all well-detected 
cases, from high values on the kinematically approaching major 
axis to low values on the receding major axis. This effect can- 
not be explained by a symmetric planar field geometry. As the 
inclination of the target galaxy increases, a pair of local maxima 
in polarized intensity begins to separate from the approaching 
major axis and propagates toward the minor axis. This is sim- 
ply a geometrical effect because maximum polarized emission 
emerges from magnetic fields perpendicular to the line of sight, 
which, in the c ase of a planar fi eld geometry, are strongest near 
the minor axis (IStil et al.ll2009t) . Even when the inclination has 
become quite substantial, and the two local maxima of polarized 
intensity are near the minor axis, there is still a systematic ten- 
dency for the receding major axis to have the overall minimum 
polarized brightness. 

This pattern in the distribution of polarized flux in our tar- 
get galaxies is visible first of all in Fig. 4 of Paper II, which 
shows maps of the polarized flux in each galaxy. We quantify 
this trend in Fig.[T]by plotting the average peak polarized inten- 
sity, < P >, and the associated Faraday depth, <<!)>, within 
inclination-corrected wedges spanning 10 degrees of azimuth in 
the disk of each galaxy in which extended polarized emission 
was detected. Azimuth is measured counter-clockwise from the 
kinematically receding major axis (see Table 1 of Paper II for 
basic data on each target). Error bars on the points in the plots 
represent the RMS variation within each wedge. The estimated 
foreground contribution to < O > and its error caused by our 
own Galaxy toward each target is indicated by the horizontal 
lines. The pairs of panels in the figure have been ordered by in- 
creasing galaxy inclination (top to bottom) extending from less 
than 10 degrees for NGC 628 to about 85 degrees in NGC 4631. 

The simple trend described above is clear, the polarized in- 
tensity showing a global minimum toward the receding major 
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axis and there being a systematic progression from one broad 
maximum at the approaching major axis to a pair of maxima that 
move to the minor axes as the inclination increases. The latter 
effect is due to geometry of a planar field, while the former one 
has no current explanation. We note that the polarization asym- 
metry along the major axis is not as pron ounced at frequencies 
of 5 GHz and higher (see e.g. lBeckl (l2007l) V 

3.2. Faraday depth 

The variation in Faraday depth with azimuth (shown in Fig. [TJ 
also shows systematic trends with increasing inclination, al- 
though not as cleanly as those seen in polarized intensity. One 
complication is that each Faraday depth distribution can be ei- 
ther positive or negative, so both options need to be considered 
in assessing a possible trend with inclination. Another compli- 
cation is the uncertainty in the foreground contribution to 
which is critical in distinguishing peaks (be they positive or neg- 
ative) in the modulation patterns from minima (consistency with 
that foreground level). What is immediately clear from an as- 
sessment of the < 0(0) > patterns is that they are in no case 
consistent with being symmetric sinusoids (either of period 2n 
or 7r) in their excursions about the estimated foreground Faraday 
depth. The implication is that a thin axisymmetric or bisymmet- 
ric spiral disk is not a viable model for the medium responsible 
for the Faraday rotat ion of the po larized emission detected at 
these frequencies (cf. lKrausSi 19901) . 

The < ^{(p) > pattern that applies to a large fraction of 
our sample has a minimum Faraday depth (consistent with the 
estimated foregound) that occurs close to the approaching ma- 
jor axis (at an azimuth of 180°) and a single maximum excur- 
sion near the receding major axis at low inclinations. This same 
pattern applies to many of the eight lowest inclination galaxies 
and would also apply to NGC 6946 if the previously published 
est imate of the Galactic foreground value, Ofc - 40radm"^ 
bv lBeckl (l200l) . were the correct one, rather than the O^c - 
23 rad m"'^ ± 2 we estimate in Pape r II. A similar considera- 
tion appUes to NGC 5194, for which iHorellou et al.l (Il992h es- 
timated <t>fc - -5radm"^ + 12 rather than our estimate of 
<t>Fc - -Hl2radm"^ + 2. We plotted these alternate estimates of 
<i>f(j in the figure with horizontal lines. We also note here a ty- 
pographical error in the value of <i>fG for NGC 4321 in Table 1 
of Paper II, which is -7radm"^ and not -17 rad m"^ as stated 
there. This pattern of maximum and minimum is not a simple si- 
nusoid but has a clear asymmetry about an azimuth of 180° that 
is most obvious in NGC 4321 and 6946. When the inclination 
exceeds 65°, there is a sudden change to a pattern of two peaks 
near the minor axis that is found in all four of the highly inclined 
galaxies in our sample. It may well be of further significance that 
several (and possibly all) of these highly inclined galaxies had 
already been identified as having a morphology suggestive of a 
circum-nuclear or galactic wind outflow. 

Another critical observation is the magnitude of the Faraday 
depth excursions from the foreground value, which is in all 
cases very modest, typically between 10 and 30 rad m"^. This 
is substantially less than the Faraday depth variations mea- 
sure d throu gh the entire disk of the Large Magellanic Cloud 
by iGaensler et al.l (i2005h using distant background sources that 
show average excursions of plus and minus 50 rad m"^ and 
peak excursions of H-245 and -215 rad nT^. Since the LMC 
is in no way a remarkable galaxy in terms of its likely field 
strength (based on the synchrotron surface brightness) or ther- 
mal electron populations (based on the star-formation-rate den- 
sity) relative to our sample galaxies, the implication is that the 



medium responsible for the Faraday rotation of the diffuse po- 
larized emission in our targets only extends over a small fraction 
of th e complete line-of-sigh t. A similar conclusion was reached 
by Berkhuiisen et alJ (119971) for NGC 5 194 based on the smaller 
Faraday depths seen at 1.4-1.8 GHz relative to 5-10 GHz for that 
target. 

An important conclusion that follows directly from a com- 
parison of the < 0(0) > and < P{(p) > plots is that the polarized 
intensity and the Faraday rotation in the highly inclined galax- 
ies of our sample must originate in distinct regions along the 
line-of-sight. This is because O is proportional to By, while P 
is proportional to B\^". For any given field geometry, a peak in 
B|| will be accompanied by a minimum in B^^" and vice versa. 
The observation that maximal excursions of both < O > and 
< P > are seen in all four of the most inclined galaxies cannot 
be achieved with any co-extensive geometry of the emitting and 
dispersing media. 



3.3. Depolarization 



In addition to detecting a systematic pattern of Faraday depth 
excursions across the full LMC disk, Gaensler et al. (2005) also 
document the very substantial depolarizing effect of the LMC 
disk on background polarized sources observed at 1.4 GHz. 
Despite a likely mean angular source size of only about 6 arc- 
sec, which projects to 1.5 pc at the LMC disk, these sources are 
more depolarized (by a factor of more than two) when the LOS 
is associated with thermal electrons in the LMC disk exceed- 
ing an emission measure of about 50 pc cm"^. The implication 
appears to be that significant RM fluctuations are present at the 
1.4 GHz observing frequency on scales << 1.5 pc. Since even the 
diffuse ionized gas that permeates galaxy disks has an emission 
measure in e xcess of about 20 pc cm ~^ in moderately face-on 
systems (e.g.. iGreenawalt et al.lll998h . we can expect a signifi- 
cant degree of depolarization for any "backside" emission in our 
galaxy sample. This will be true in particular for the diffuse po- 
larized emission originating in the target galaxy itself, since the 
relevant scale is then that of the observing beam, >15 arcsec, 
which projects across more than 700 pc at the typical galaxy dis- 
tance of 10 Mpc. 

It therefore seems likely that the polarized intensity from low 
inclination galaxy disks observed at GHz frequencies is domi- 
nated by emission from that portion of the galaxy disk/halo that 
faces us. Corresponding structures from the far side of the galaxy 
would be dispersed and depolarized by turbulent magneto-ionic 
structures in the star-forming mid-plane. We recall that it is only 
in those regions dominated by a regular rather than a turbulent 
magnetic field that a net polarized emissivity is expected at all. 
Given the strong concentration of massive-star formation and its 
associated turbulent energy injection into the mid-plane, there 
may well be two zones of enhanced net polarized emissivity that 
are offset above and below the mid-plane. Each region of polar- 
ized emissivity will then experience the dispersive effects of its 
own line-of-sight foreground, which is likely to be dominated by 
thermal rather than relativistic plasma. For the near-side polar- 
ized emission, this is likely to be caused by the extended thermal 
halo of the host galaxy, while for the far-side polarized emission 
there will be the additional contribution of the dense (depolariz- 
ing and dispersive) mid-plane. 
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3.4. The "second" polarized disk 

Because of the likely Faraday depth of the mid-plane medium, 
it is conceivable that the near- and far-side polarized emission 
zones in moderately face-on galaxies would experience very dif- 
ferent Faraday dispersion, making it possible to distinguish the 
two components along each line-of-sight usingthe RM synthe- 
sis technique, which we employed in our study. We conducted 
a deep search for multiple Faraday depth components along the 
line-of-sight to detect disk emission from each of our sample 
galaxies, after applying a spatial smoothing to the P{<p) cubes 
that results in an angular beamsize of 90 arcsec. Solid detec- 
tions of faint secondary emission components were made in the 
brightest face-on galaxies of our sample, as shown in Figs. |2]- 
|4] In addition to the bright polarized emission that typically re- 
sides at a Faraday depth within only a few tens of rad m"^ of 
the Galactic foreground value, much fainter polarized emission 
(by a factor of 4 - 5) is detected in NGC 628, 5194, and 6946, 
which is offset to both positive and negative Faraday depths by 
about 200 rad m"^. We first considered whether these faint sec- 
ondary components might be instrumental in nature, since they 
are similar to the Faraday depth side-lobes of our instrumental 
response (as described in Paper II), but concluded that they are 
likely reliable. Indicators of the reliability of these faint com- 
ponents are that (1) they do not occur toward the brightest low 
dispersion components; (2) the Faraday depth separation of the 
secondary components varies from source to source, while the 
instrumental sidelobe response does not; and (3) the faint posi- 
tive and negative -shifted components form a complementary dis- 
tribution of eachother, rather than merely repeating themselves 
in detail. 

Detection of this highly dispersed and probably depolarized 
"second disk" supports the emerging model in which the polar- 
ized intensity observed at GHz frequencies from nearby galaxy 
disks is dominated by a region of emissivity located above 
the mid-plane on the near-side, which subsequently experiences 
Faraday rotation within the extended halo of the galaxy amount- 
ing to only a few tens of rad m"^. Significantly fainter polarized 
emission (by a factor of 4 - 5) is detected from the farside of 
the mid-plane, which displays the additional dispersive effects 
of that mid-plane zone amounting to plus and minus 150 - 200 
rad m"^ within the three galaxies where this could be detected. 



3.5. Small-scale RM fluctuations 

Extremely small-scale RM fluctuations within our own 
Galaxy were discovered in our sample data for the field of 
NGC 7331. This target is observed through the Galaxy near 
{l,b) ~ (94,-21), within an extended (30° X 30°) region 
of particularly large negat ive RMs near -200radm"^ (e.g., 
iJohnston-Holhtt et al.l 12004}) . A corresponding region of large 
positive Galactic RMs is centered near (l,b) ~ (250,-10). 
These two regions correspond to the directions where we look 
directly along what is likely to be t he axisymmetric spiral field 
of the Galaxy (e.g.. lSun et alJl2008l) . Diffuse polarized emission 
from the lobes of a background head-tail radio galaxy, the disk 
of NGC 7331, and even the Galactic synchrotron itself display 
the remarkable oscillatory behavior of the Faraday depth of the 
polarized emission with position in the field. As shown in the 
Faraday depth versus declination slices of Figs. |5] and |6] there 
are two dominant Faraday depths present in this field, one near 
-ISOradm"^ and the other near Oradm"^. Depending on the 
exact location along the indicated declination slice, either one 
or the other of these RMs are encountered. In some regions, the 



transition from one RM to the other is well-resolved and a sin- 
gle RM value is observed over several beam-widths (of about 
20 arcsec or 0. 1 pc at a distance of 1 kpc), while in other re- 
gions the transition is unresolved, such that both RMs overlap 
spatially with only the peak polarization showing the oscilla- 
tion between the two values. When completely unresolved extra- 
galactic sources are observed in this field, they display only one 
or the other of these two possible foreground RMs (compare 
Table 2 of Paper II), but more extended sources show an oscilla- 
tion between the two values, - 1 80 rad m"^ occurring 3-4 times 
as often as Oradm"^. The regular Galactic magnetic field ap- 
pears to be organized into discrete filamentary components with 
transverse sizes significantly less than a parsec. 

4. Magnetic field distributions 

The observational results summarized in Sect.[3] suggest a mag- 
netic field geometry dominated by an in-plane field described 
by an axisymmetric spiral structure (ASS) and/or possibly a bi- 
symmetric spiral structure (BSS). These topologies assume a 
planar field configuration (in the ix,y) plane) given by a fam- 
ily of logarithmic spirals defined in radial coordinates (r, (f>) by 
the form 



M+c) 



(2) 



for a radial scaling constant, a, an angular scaling constant, b, 
which is related to the spiral pitch angle, t//[^„ by. 



-a;, = tan-'(/7) 



(3) 



and an angular offset, c, which defines each curve in the family. 
We find it convenient to use the complement of the spiral pitch 
angle, i^ i,, given by. 



= 90°-f;,, = tan-'(l/^7). 



(4) 



The right-handed cartesian components of an inward directed 
planar ASS logarithmic spiral are then given by. 



B.v = B cos(0 -I- i/r.vy) 
By = Bsm(4) + il/^y) 
B- = 0. 



(5) 
(6) 
(7) 



For the corresponding case of a planar BSS magnetic field with 
the usually assumed sinusoidal modulation of B with 0, we have. 



Bx - B COS(20 + ij/xy + ^A 

By = B sin(20 -H i/r^v + /") 
B, = 0, 



(8) 
(9) 
(10) 



where fi is used to track the positive peak in the field modulation 
pattern from an initial value fiQ such that. 



fi^luo + [ln(r/a)]/b 



(11) 



at radius r, in terms of the spiral scaling constants defined above, 
a and b. 

The out-of-plane field topology is not well-constrained by 
previous observations, but might be expected to have either an 
even or odd configuration o f the symmetry about the mid-plane 
(see Fig. 2 of lWidrowll2002l) . An even configuration corresponds 
to the case where the azimuthal (toroidal) component of the field 
has the same sign both above and below the mid-plane. The re- 
sulting field geometry has a quadrupole structure in the poloidal 
field and is the one predicted to emerge most naturally from an 
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ao) dynamo operating at intermediate to large rad ii in a galac- 
tic di sk, where differential rotation is important (Elstn er et al.l 
II992I) . An odd-parity configuration has the opposite signs of 
the azimuthal field above and below the mid-plane. The asso- 
ciated structure of the poloidal field is that of a dipole. This 
topology may be associated with the dynamo process that 
may be dominant within the circum-nuclear regions where solid 
body rotation may prevail in galaxies (El stner et al. 1992). The 
aoj dynamo in galactic halo s may also generate dipolar fields 
ISokoloff & Shukurovl (Il990h . 

Of stronger relevance is the expected height above and below 
the mid-plane at which the polarized synchrotron emission might 
originate. From the discussion in Sect.[3] we expect that polar- 
ized emissivity may peak on either side of the galaxy mid-plane, 
but that the near-side component will dominate the detected po- 
larized intensity at GHz frequencies because of depolarization 
in the turbulent mid-plane, which is an intervenor to the far-side 
emission. The Faraday depth distribution will reflect all of the 
relevant propagation effects that the emission has experienced. 
For the near-side component this will reflect the extended near- 
side halo of the target galaxy, while for the far-side component, 
the additional dispersion in the mid-plane region wiU also con- 
tribute. 

We extend the usual planar ASS and BSS field topology 
with the addition of an out-of-plane component taken from a lin- 
ear combination of dipole and quadrupole topologies. The basic 
equations describing dipole a nd quadrupole field s and magnetic 
flux functions can be found in lLong et al.l (12007!') . For the simple 
case under consideration, there is cylindrical symmetry about 
the galaxy rotation axis, z. In terms of the angle from the rota- 
tion axis, 0, and the distance from the origin, p, the two perpen- 
dicular components of the poloidal magnetic field with a dipole 
moment, D, and quadrupole moment, Q, (each with a non-trivial 
sign) are given by. 



Bp = 



2Dcos{0) 3Q[3cos\6) - 1] 



Dsin(6) 3Qsin(6)cos(6) 

Be - 5 — + T-j ■ 

pi 2p4 

The corresponding magnetic flux function, M^, is given by, 
p Ap^ 



(12) 
(13) 

(14) 



The surfaces defined by 'I' = constant represents example sur- 
faces on which the magnetic field lines reside. The field compo- 
nents out of- and within the plane are given by. 



B, = BpCos(_e)-Bgsin(e), 
Br = Bpsin(e) + Bgcos(6), 



(15) 
(16) 



which yield the total field strength and local orientation angle 
from. 



B^ = B^ + B? 



i/r, = tan"' (B,/Br). 



(17) 
(18) 



The poloidal field topologies are shown in a plane that includes 
the rotation axis in Fig. |7] for pure dipole, quadrupole, and a 
mixed dipole plus quadrupole configuration. The units of the two 
axes in the plot are arbitrary, since the topologies are self-similar. 
A three dimensional depiction of the pure dipole and quadrupole 
field topologies is given in Fig. [8] where several surfaces of con- 
stant magnetic flux are shown for each case. 



Combining the planar (Eq. |4]i and out-of-plane (Eq. HM ge- 
ometries permits the three right-handed cartesian components of 
the modified ASS magnetic field in the frame of the galaxy to be 
written as. 



B f = B cos(0 + i/'vy) cos(i^^). 
By = B sin(0 + ifr„) cos(i/'-;), 
B, = Bsin(i/r,). 



(19) 
(20) 
(21) 



When viewed at an inclination, /, this yields the observer's frame 
components, 

Bf' - B cos(0 + t//„,) cos(i^,), (22) 
By' - B sin(0 + i^fv) cos(i^,) cos(/) - B sin(i/'j.) sin(/), (23) 
B,' - B sin(0 + i^.ty) cos(i^,) sin(0 + B sin(i^,) cos(0, (24) 

where the (x',y',z') are the major axis, minor axis, and line 
of sight, respectively. If the spiral is a trailing one (as demon- 
strated in every kinematically studied galaxy with a well-defined 
nearside) then the positive x' axis so defined corresponds to the 
receding major axis. The projected parallel and perpendicular 
components of the magnetic field and the orientation angle of 
Bx are given by. 



Bii 



B-, 



+ Bl 



Xq - arctan 



/B, 



(25) 
(26) 

(27) 



We recall that the Faraday depth caused by a magneto-ionic 
medium is proportional to By, the polarized intensity is propor- 
tional to B^^", and the intrinsic polarization angle (giving the 
E-field orientation) i^xvi - x'q~ 90°. 

In the corresponding case of a modified BSS magnetic field, 
we have 



Bx - B cos(20 -I- ijjxy + i-i) cos((/'j.). 
By = B sin(20 + + /j) cosiiff^), 
B, = Bsin(i/r,). 



(28) 
(29) 
(30) 



When viewed at an inclination, this yields the observer's frame 
components, 

B y, = B cos(20 -I- if/„. + fi) cos(i/';), (31) 
By' - B sin(20 + ^„ + lA cos(tff,) cos(/) - B sin{ilf,) sin(0, (32) 
B-, - B sin(20 + i/Zxy + n) cos(i^,) sin(/) -1- B sm(if/^) cos(/), (33) 

and the corresponding By, B^, and;\fg as above. 

With these definitions in place, it is possible to explore the 
parameter space of modified ASS and BSS field topologies and 
produce both images and azimuthal traces of the expected distri- 
butions of B|| and B^. We use these measures as proxies of the 
Faraday depth and polarized intensity, respectively. In the ab- 
sence of a model for the spatial distributions of cosmic-ray and 
thermal electrons, we do not attempt to reproduce the observ- 
ables directly, but merely the systematic modulations of P and O 
with azmiuth. 

We present images and traces of the ASS and BSS predic- 
tions in Figs. l9l-[T2l Each group of 2 x 2 panels in the figure 
represents a contour plot of By and B^ (at top) and a series of 
traces of B||(0) and B^{(p) at fixed radii (below). The relative 
weights of a dipole (D) and quadrupole (Q) field (each includ- 
ing a nontrivial sign), galaxy inclination, /, spiral pitch angle. 
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ipxy (positive for counter-clockwise and negative for clockwise), 
and distance(s), Z or (Zl and Z2), from the mid-plane (positive 
toward the observer and negative away) are indicated at the top 
of the contour plots, together with an ASS or BSS designation. 
Two representative spirals are drawn for reference, both directed 
inwardly for the ASS plots and one inward and the other out- 
ward for the BSS plots. The inclination is defined such that for 
kinematically trailing spirals, the receding and approaching ends 
of the major axis are as indicated. The "near" and "far" designa- 
tions in the plots also aid in demonstrating the spatial orientation 
of the disk unambiguously. An ellipse that is offset from the disk 
by the distance Z is drawn to demonstrate its location relative to 
the mid-plane. The radii at which the azimuth plots were made 
are marked by the same linetype on the contour plots. The az- 
imuth angle increases counter-clockwise from the receding ma- 
jor axis. 

We also plot the projected angle of B^, x[M')^ within the 
same panel that presents the trace of B^{(p) using the right-hand 
scale. ThexQW traces are only plotted for the ASS cases, since 
the BSS cases vary so dramatically with radius. 

For simplicity, we begin with a fixed spiral pitch angle, 
(/''j,, - 20°, since this corresponds we ll with the avera ge mea- 
sured value for our sample galaxies (cf. lKennicutlll981h . which 
vary from about 15-25°, and we begin by considering counter- 
clockwise (CCW) spirals. Other values will also be considered 
below. First we consider the predictions for the simple planar 
ASS and BSS spirals in Fig|9] A radial spiral scaling constant, 
a = 5, and a galaxy disk radius, r,„a j - 35, are assumed for 
illustration. These choices have no effect on the results. 

4.1. Planar models 

Planar, axisymmetric field topologies yield projected field com- 
ponents (as shown in the left-hand groups of Fig. ^ that are 
independent of radius and have very simple symmetries, B||(0) 
having positive and negative excursions that are fully symmetric 
about zero. The single positive peak occurs at ~ ifr'^^. offset 
from the approaching major axis, the negative peak being oppo- 
site to this (as also demonstrated by Krause 1990). The Bj.(0) 
component has a complimentary behavior with two equal min- 
ima slightly offset from both major axes, and two equal maxima 
slightly offset from both minor axes. As the inclination increases 
from face-on toward edge-on, the amplitude of the B||(0) and 
B±{(p) modulation increases, although there is no change in the 
azimuthal location of either the maxima or minima. The pro- 
jected orientation of the plane-of-sky field, XoW^ has the ex- 
pected linear variation with azimuth for a nearly face-on geom- 
etry, which becomes increasingly non-linear as the inclination 
increases. 

Changing the sense of the field to be outwardly directed, 
rather than inwardly directed, changes only the sign of B||(0), 
and leaves Bj.(0) unchanged. This is true for all of the field ge- 
ometries we consider 

The basic BSS projected field patterns and their variation 
with galaxy inclination are illustrated in the right-hand groups 
of Fig. |9l Because of the modulation of field sense with azimuth 
along the family of spirals, there is a strong radial dependence 
of the projected field components and their variation with az- 
imuth. B\\((f)) has positive and negative excursions that are still 
fully symmetric about zero, but it exhibits two equal maxima 
and two equal minima (rather than only one for ASS), which 
are located at different azimuth depending on the radius. Bj.(0) 
displays four equal maxima at the azimuthal angles where B||(0) 
has its zero crossings, and four equal minima in-between (at the 



locations of maximum and minimum By). As with the basic ASS 
fields, a higher inclination yields an increased fractional modula- 
tion of B||(0) and B±{(p), while not affecting the location of those 
excursions. The strong radial dependence of these azimuthal pat- 
terns implies that if there were significant averaging of radii (by 
even as little as 20% in Ar/r for a pitch angle of i^'j.,, ~ 20°) 
then most of the predicted modulation would disappear This is 
particularly true of Bj.(0), for which the predicted modulation is 
both intrinsically smaller and twice as rapid as that of B\\{(p). The 
strong radial dependance of BSS models severely limits their 
predictive power for the general patterns of azimuthal modula- 
tion that are observed in galaxies and so will not be considered 
further 

The azimuthal modulation of By and B± in the planar ASS 
and BSS spiral topologies is apparently inadequate for describ- 
ing the general patterns noted above in Sect.[3] Neither of these 
topologies predict a clear distinction between the receding and 
approaching kinematic major axis in terms of either the intrinsic 
brightness of polarized intensity or the Faraday depth that it may 
encounter 

4.2. Thick disk models 

In Figs.[T0l-fT2l we explore the extension of the basic ASS topol- 
ogy with either a dipole (D = 1, left-hand groups) or quadrupole 
(Q = 1, right-hand groups) modification to the out-of-plane field, 
as given by Eq.|2T| In keeping with the discussion of Sect.[3] we 
calculate a projected distribution that is populated within a spec- 
ified range of mid-plane heights extending from Zl to Z2 on the 
near-side of the galaxy and from -Z2 to -Zl on the far-side. For 
lines-of-sight that intersect the mid-plane disk within the nomi- 
nal galaxy radius, only the near-side region (Z > 0) contributes to 
the integral, to approximate the effect of mid-plane depolariza- 
tion. For other lines-of-sight, out to 1 .2 times the nominal galaxy 
radius, both the nearside and farside zones are integrated. The in- 
tegral is determined by first evaluating By and B^ at a sequence 
of finely-sampled line-of-sight depths and then forming the aver- 
age. The two lower panels of each group were modified to allow 
more direct comparison with the data of Fig. [Uby plotting the 
mean modeled quantity within azimuthal wedges rather than the 
azimuthal traces at discrete radii. The error bars represent the 
RMS variation within each wedge, which are principally caused 
by a dependence on radius. 

We first consider the dipole and quadrupole models that are 
integrated over mid-plane heights of Izl - 2 — > 4 (about 10% 
of the disk radius) in Fig. [TO] Although the azimuthal trace of 
projected field modulation with azimuth show some variation in 
amplitude with radius (as indicated by the error bars), there are 
well-defined maxima and minima. The local modification to the 
planar spiral field is illustrated in Fig. Q by the vector orienta- 
tions along the horizontal line drawn at Z = +4. The addition of 
a poloidal field component has introduced several notable differ- 
ences from the planar case. Although B||(0) exhibits maxima and 
minima at the same azimuth as previously (near the major axes), 
it no longer has excursions that are symmetric about By - 0, 
By(0) being offset to negative values of By (for an inwardly- 
directed spiral field) and showing complementary modulation in 
the dipole and quadrupole cases. The character of Bj^(0) is more 
substantially modified by the addition of a poloidal field. We first 
consider the quadrupole case (Q = 1) shown in the right-hand 
groups of Fig.fTOl 

At low inclinations, a single minimum in B^ occurs near the 
receding major axis at ~ t//\^, and a broad single peak centered 
near the approaching major axis. As the inclination increases. 
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the single global minimum in Bj^ remains close to the receding 
major axis, while the single broad maximum divides into two 
maxima that separate and shift toward the two ends of the mi- 
nor axis. This is exactly the trend shown by the data in Fig. [1] 
The other noteworthy attribute of these distributions is that the 
largest excursion in |B||| is found on the receding major axis at 
all inclinations, representing the greatest Faraday depth at this 
azimuth for coextensive emitting and dispersive media. For low 
inclinations, there is a broad minimum in |B||| at the approaching 
major axis, which becomes a secondary peak at this location (of 
opposite sign) as the inclination increases. This roughly agrees 
with the patterns seen for the majority of lower inclination galax- 
ies in Fig. [1] these models still being more symmetric than the 
data, which exhibit a much clearer nearside/farside asymmetry. 

For the dipole case (D = 1), shown in the left-hand groups of 
Fig. [To] the receding and approaching major axes are reversed, 
since the poloidal field component of the dipole at positive mid- 
plane offsets is directed outward and not inward (as shown in 
Fig. |7]l. For this case, it is the approaching major axis that is 
predicted to have a global minimum in B^, and hence also the 
minimum intrinsic polarized intensity. In a similar way, it is the 
approaching major axis for which the greatest value of 
hence the largest associated Faraday depth is found. The degree 
of modulation with azimuth is rather modest for the dipole case 
integrated over these heights. 

From this comparison, it is clear that a large-scale 
quadrupole ASS field topology provides an excellent model for 
explaining the modulation of polarized intensity with azimuth, 
while the dipole clearly does not. The modulation of Faraday 
depth with azimuth is also explained qualitatively for the low 
inclination systems, but is not reproduced in detail. The crucial 
attributes of the quadrupole model in matching the observed az- 
imuthal modulation pattern, < P(0) >, are the magnitude and 
particularly the sign of of the angle tfr^ as given by Eq. [18] For 
typical values of the spiral pitch angle, \if/'^^.\ < 25°, the projected 
field component. By of Eq. [24] has minima and maxima near 
= and 180°. The "hour-glass" shape of the quadrupole field 
yields a positive sign of i//, for small positive values of z, i.e. 
toward the observer, and this is what yields a minimum of B± 
near = for the nearside emission. The "donut" shape of the 
dipole field, on the other hand, yields a negative sign of 0, for 
small positive values of z, resulting in a maximum of B^ near 
= 0. 

4.3. Halo models 

We now consider dipole and quadrupole models that extend to 
more substantial mid-plane heights in an attempt to reproduce 
the azimuthal modulations of Faraday depth seen in the data. 
In Fig. [TT] we illustrate the result of integrating the same mod- 
els just considered over mid-plane heights of Izl = 4^10 
(or about 30% of the disk radius). By considering first the 
quadrupole models in the right-hand panels, we see substantial 
similarity between these models and those originating closer to 
the mid-plane (Fig. [TOl l. The modulation pattern in Bj.(0) has a 
slightly smaller amplitude. The modulation pattern of the halo 
B||(0), while similar to its thick disk counterpart, has a clear 
asymmetry in the approaching major axis minimum about an 
azimuth of 180°, which is reminiscent of what is seen in the data 
of Fig. [T] The pattern shown is for a positive pitch angle of 20° 
(CCW spiral), while for a negative pitch angle (CW spiral) it is 
mirrored about an azimuth of 180°. We recall that changing the 
sign of the spiral field from inward directed (as shown) to out- 
ward directed changes only the sign of B||(0). Qualitative agree- 



ment of the halo predictions for the quadrupole with the basic 
< <I>(0) > patterns is reasonable for many of the low inclination 
galaxies of our sample. The dipole halo models (shown in the 
left-hand panels) with their maximal excursion near an azimuth 
of 180° reproduces the observations less successfully. 

Although the models considered to this point can reproduce 
the general pattern of < P(0) > modulation at all galaxy incli- 
nations and the < '^(<p) > modulations for many low inclination 
targets, they have not provided any agreement with the distinc- 
tive doubly peaked pattern of < O(0) > around the minor axes 
seen in our four highest inclination galaxies. We now consider 
an additional variant of our halo models in an attempt to repro- 
duce those < '5(0) > patterns. In Fig.[T2]we consider dipole and 
quadrupole models for an extremely large pitch angle of 85°, to 
illustrate the impact of a planar field component that is essen- 
tially radial. The distinctive doubly-peaked pattern of < O(0) > 
at large inclinations can be reasonably reproduced in this case, 
but only for the dipole topology. 

5. Discussion 

Our analysis of the projected three-dimensional magnetic field 
topologies presented in Sect. |4] and their predicted observable 
consequences for the azimuthal modulation, B||(0) and B^(<p), 
has provided a plausible explanation of the very general ob- 
served trends noted in Sect. [3] A self-consistent scenario has 
emerged that accounts for the polarized intensity and its Faraday 
dispersion observed at GHz frequencies from galaxy disks. The 
detected polarized intensity is dominated by a zone of emissivity 
above the mid-plane on the side of the galaxy facing the observer 
(at a height of perhaps 5 - 10 % of the disk radius). This thick 
disk emission arises in a region that is dominated by an axisym- 
metric spiral with an out-of-disk quadrupole topology, which is 
responsible for a distinctive modulation of B^{(p) and its varia- 
tion with galaxy inclination. This emission is affected by only a 
modest amount of Faraday dispersion, of a few tens of rad m"^, 
within the nearside halo of the galaxy in its subsequent prop- 
agation. For the majority of low to modest inclination galaxies 
(< 60°), the dispersive foreground topology is consistent with an 
extension of the thick disk ASS quadrupole out to larger heights 
above the mid-plane (of perhaps 30% of the disk radius). 

The most highly inclined galaxies of our sample require 
an alternative halo field topology, in the form of a radially- 
dominated dipole, which yields a distinctive doubly-peaked 
modulation of '5(0). It seems significant that in many or pos- 
sibly all of the highly inclined galaxies of our sample there 
is evidence of a significant circum-nuclear outflow component 
to the polarized emission, in addition to that of the disk. This 
circum-nuclear component would quite naturally be expected to 
be associated with a dipole, rather than a quadrupole field, in 
view of the likely dominance of the over the a cL) dynamo 
mechanism at small galactic radii (e.g. lElstner et al.i ri992). This 
circum-nuclear dipole field would also be less likely to have any 
association with the spiral pitch angle of the disk given its ori- 
gin. Because of the shallower roll-off' with radius of a dipole 
compared to the quadrupole field (see Eq. [TSJ, a dipole com- 
ponent may come to dominate the h alo field of the associated 
galaxy when both are present. ISokoloff & Shukurovl(ll990l ) also 
argued that a ato dynamo operating in the halo would directly 
produce a dipole field. Non-stationary global halo models (e.g. 
Brandenburg et al.'l992) may also provide a natural explanation 
of the dipole signature on the largest scales. 

In addition to the bright polarized emission originating in 
the nearside, the corresponding rear-facing region of polarized 
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emissivity of the thick disk can also be detected in relatively 
face-on galaxies if sufficient sensitivity is available. This emis- 
sion is substantially weaker, by a factor of 4 - 5, and consistent 
with depolarization caused by fluctuations in the magneto-ionic 
medium of the mid-plane on scales smaller than a pc. This fainter 
polarized component is affected by much greater Faraday disper- 
sion, corresponding to both plus and minus 150 - 200 rad m"^ 
in its propagation through the dense mid-plane plasma, as well 
as the near-side halo. These two maxima (one positive and one 
negative) in Faraday depth are aligned approximately with the 
major axes of each galaxy, and have approximately symmetric 
excursions about the Galactic foreground value. This pattern is 
consistent with the expectation for a simple planar ASS field in 
the galaxy disk. 

Future observations of nearby galaxy disks at frequencies 
below 200 MHz, such as with the upcoming LOFAR facility, 
will likely detect net polarized emission from even larger heights 
above the galaxy mid-plane and exclusively from regions unob- 
structed by the mid-plane in projection. A good indication for 
the predicted observables is given in Fig.[TT]in which we present 
model integrations of the upper halo (out to 30% of the disk ra- 
dius). 
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Fig. 1. Azimuthal variation in peak polarized intensity (center panel) and associated Faraday depth (left panel) for galaxies with 
extended polarized emission (illustrated in the right panel from Fig. 5 of Paper II). The mean values in azimuthal wedges, each 
subtending 1 0° within the galaxy disk, are plotted with error bars giving the wedge RMS . Galaxies are arranged in order of increasing 
inclination (top to bottom) from face-on. Azimuth is measured counter-clockwise from the receding major axis. 
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Fig. 1. (continued) Azimuthal variation in peak polarized intensity (center panel) and associated Faraday depth (left panel) for 
galaxies with extended polarized emission (illustrated in right panel). 
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Fig. 1. (continued) Azimuthal variation in peak polarized intensity (center panel) and associated Faraday depth (left panel) for 
galaxies with extended polarized emission (illustrated in right panel). 
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Fig. 2. Polarized intensity at distinct Faraday depths toward NGC 628. The dominant Faraday depth component, centered near 
-30rad m~^ is shown in the center panel, while the two secondary components centered near -213 and +145 rad m~^ are shown on 
the left and right. The greyscale varies as indicated. The contours begin at 0.29 mJy beam"^ and increase by factors of 1.1 for the 
secondary components and 1.3 for the primary component. 
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Fig. 3. Polarized intensity at distinct Faraday depths toward NGC 5194. The dominant Faraday depth component, centered near 
+ 13 rad m"^ is shown in the center panel, while the two secondary components centered near - 1 80 and +200 rad m"'^ are shown on 
the left and right. The greyscale varies as indicated. The contours begin at 0.6 mJy beam"^ and increase by factors of 1.1 for the 
secondary components and 1.3 for the primary component. 
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Fig. 4. Polarized intensity at distinct Faraday depths toward NGC 6946. The dominant Faraday depth component, centered near 
+38 rad m"^ is shown in the center panel, while the two secondary components centered near -162 and +228 rad are shown on 
the left and right. The greyscale varies as indicated. The contours begin at 0.7 mJy beam ' and increase by factors of 1.1 for the 
secondary components and 1.3 for the primary component. 
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Fig. 5. Polarized intensity as a function of Faraday depth and declination in the field of NGC 733 1 . The two left hand panels are at 
right ascensions that intersect the diffuse lobe of a background head-tail radio galaxy. The right-hand panel intersects a background 
double-lobed radio galaxy as well as a region of diffuse Galactic polarized emission. 
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Fig. 6. Polarized intensity as a function of Faraday depth and declination in the field of NGC 7331. Slices are presented at three 
right ascensions that intersect the disk of NGC 7331 near the major axis. 
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Fig. 7. Depiction of the assumed poloidal modification to the out-of -plane field topology. A planar logarithmic spiral is modified 
by the local orientation of a dipole or quadrupole field that is symmetric about the rotation axis. The panels are labelled with the 
relative strengths of dipole (D) and quadrupole (Q) moments and illustrate a pure dipole, pure quadrupole, and a 1 : 100 mix of dipole 
and quadrupole from left to right. Dashed horizontal lines illustrate the heights for which model distributions are shown in Fig.fTOl 




Fig. 8. Depiction of the magnetic flux function for a pure dipole (left) and quadrupole (right) field. In our modeling we assume that 
a planar logarithmic spiral field, as illustrated with the solid lines in the Figure defines the angle 4i„, while the local orientation of 
a dipole or quadrupole field that is symmetric about the rotation axis defines the angle, ij/^, out of the plane. 
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Fig. 9. Simulated field distributions. Each group of 2 x 2 panels presents contour plots (top) and azimuthal traces (bottom) of B\\ 
(left) and B± (right). The dipole and quadrupole moments (D, Q), galaxy inclination (1), spiral pitch angle (P), and height above 
the galaxy mid-plane (Z) are indicated at the top of each contour plot. Contours are drawn at +10, +20, ±30, . . . +90% of \B\. Two 
reference spirals are drawn, both with field directed inward {B - +1) for the axisymmetric (ASS) case and one in and the other out 
for the bisymmetric (BSS) case. The approaching and receding major axis of the galaxy are indicated for trailing spirals. Azimuth 
is measured CCW from the receding major axis. Azimuth traces are drawn for the radii indicated in the contour plots. A trace of 
the orientation of B±, XqW' given in the B^ panel for ASS models with a dashed linetype using the right-hand scale in degrees. 
Here we compare the ASS (left hand groups) and BSS (right hand groups) planar models (Z - 0). The galaxy inclination is 20° in 
the upper series of groups and 40° in the lower series. The spiral pitch angle is 20°. 
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Fig. 9. (continued) Simulated field distributions. Here we compare ASS (left hand groups) and BSS (right hand groups) planar 
models (Z = 0). The galaxy inclination is 60° in the upper series of groups and 80° in the lower series. 
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Fig. 10. Simulated field distributions (as in Fig.|9ll. Here we compare dipole- (left hand groups) and quadrupole- (right hand groups) 
ASS models integrated over a zone from Zl to Z2 (|Z| - 2-4). For those lines-of-sight that intersect the mid-plane of the disk, only 
positive values of Z are included in the integral. The galaxy inclination is 20° in the upper series of groups and 40° in the lower 
series. 
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Fig. 10. (continued) Simulated field distributions. The galaxy inclination is 60° in the upper series of groups and 80° in the lower 
series. 
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Fig. 11. Simulated field distributions (as in Fig.|9ll. Here we compare dipole- (left hand groups) and quadrupole- (right hand groups) 
ASS models integrated over a zone from Zl to Z2 (|Z| - 4-10). For those lines-of-sight that intersect the mid-plane of the disk, only 
positive values of Z are included in the integral. The galaxy inclination is 20° in the upper series of groups and 40° in the lower 
series. 
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Fig. 11. (continued) Simulated field distributions. The galaxy inclination is 60° in the upper series of groups and 80° in the lower 
series. 
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Fig. 12. Simulated field distributions (as in Fig.|9ll. Here we compare dipole- (left hand groups) and quadrupole- (right hand groups) 
ASS models integrated over a zone from Zl to Z2 (|Z| = 4-10) with a large pitch angle (P = 85°). For those lines-of-sight that 
intersect the mid-plane of the disk, only positive values of Z are included in the integral. The galaxy inclination is 20° in the upper 
series of groups and 40° in the lower series. 
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Fig. 12. (continued) Simulated field distributions. The galaxy inclination is 60° in the upper series of groups and 80° in the lower 
series. 



